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Conditions for intra and inter layer Coulomb interactions in multi-walled carbon nanotubes are 
derived from stability of capacitance excitations. It is pointed out, supposing the stability conditions 
are not satisfied, that the system has unstable modes which correspond to a charge transfer between 
layers or charge density oscillation in each layer depending on zero mode or non-zero modes. It is 
argued that the stability conditions can be broken when the vacuum polarization processes due to 
the massive bands are taken into account. 
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I. INTRODUCTION 

Electrical properties of carbon nanotubes(CNTs ||, 
^) have attracted much attention from various points 
of view. They are characterized by the quantum me- 
chanical behavior of the 7r-electrons which are interact- 
ing with each other via the Coulomb interaction. Sev- 
eral aspects of the Coulomb interaction are already seen 
in some experiments of charging and of temperature 
dependent resistivity j|] which are consistent with cur- 
rent understanding of low energy excitations in metallic 
CNTs. The low energy excitations in the metallic CNTs 
are theoretically modeled by massless bands (linear dis- 
persion bands) with the Coulomb interaction [|| ^, 
On the other hand, the massive bands (other sub-banc 
except the massless bands) modify the Coulomb interac- 
tion between massless fermions (electrons in the massless 
bands) through the vacuum polarization effect j|] . As for 
single-wall CNTs with small diameter, the correction to 
the long wave length modes of the Coulomb interaction is 
negligibly small. However in case of multi-walled CNTs, 
because of their large diameter, the gap energy of the 
massive bands is small; besides, the number of the mas- 
sive bands is large so that the correction would be large 
also. Because of that the low energy spectrum of the 
multi-walled CNTs is thought to be under the influence 
of the vacuum polarization due to the massive bands. 

The multi-walled carbon nanotubes (MWCNTs) con- 
sists of several separate single-wall CNTs and may have 
an unique electrical property which is not shown in each 
individual single-wall CNT. In this paper, we investi- 
gate a novel electrical property of a MWCNT which 
has two layers. Even in this simple case, there is ac- 
tually a different aspect as compared with the single- 
wall CNT. That is the Coulomb interaction between the 
layers: interlayer Coulomb interaction. The Coulomb 
interactions in this system consist of three potentials; 
intra-layer Coulomb interactions in each layer and the 
interlayer interaction. Although it is difficult to obtain 
information on the Coulomb potentials as compared with 
the one in the single-wall CNTs, it can be thought that 
the intra-layer Coulomb potential for an inner layer is 
the same potential in the single-wall CNTs, however the 
intra-layer one for an outer layer is not simple because 
of the screening due to the presence of the electrons in 



the inner layer. The interlayer unscreened Coulomb po- 
tential between the layers in the MWCNTs is calculated 
theoretically |]l0| . 

Since each layer of the MWCNTs has it's own heli- 
cal structure, we can divide the MWCNTs into three 
categories depending on the structural property of 
each layer. That are (1) metallic-metallic(M-M) , (2) 
metallic-semiconducting(M-S) and (3) semiconducting- 
semiconducting (S-S) MWCNT. If there were no inter- 
layer interaction in the MWCNTs, analysis of those sys- 
tem is the same as that of the two individual single-wall 
CNTs and the stability of the system is trivial. But there 
is the interlayer Coulomb interaction. This interaction 
cause a new excitation mode; charge transfer between 
the two layers: capacitance mode which may cause an 
instability of the system. For S-S and S-M types, the 
stability of the system is thought to be rigid. However 
we need to consider M-M type MWCNTs carefully. The 
Coulomb potentials decide the energy spectrum of the ca- 
pacitance mode. Hence, for some kinds of potential, the 
mode could be unstable, that is, it's excitation spectrum 
have a negative eigenvalue. This cause the instability of 
the vacuum of the system and will change the vacuum 
into an another new vacuum state that we do not know 
within the model used in this paper. It is the purpose 
of this article to find conditions on the stability of the 
system and check if the condition is satisfied including 
the vacuum polarization correction due to the massive 
bands. 

The paper is composed as follows. A model Hamil- 
tonian for the low energy excitations in the M-M type 
MWCNTs is constructed in Sec. II. In Sec. Int we study 



the zero mode spectrum of the system and establish a 
condition for the stability of the zero mode. Then we in- 
vestigate the non zero modes in Sec. |^ and check if the 
system is stable in realistic Coulomb interactions. Sum- 
mary and discussion are given in Sec. 0. 



II. LOW ENERGY EFFECTIVE MODEL 

Because the kinetic terms and the Coulomb interac- 
tions consist of densities of each layer, we begin by spec- 
ifying the densities in both layers. Let us define the 
density for the outer layer as Jout and for the inner 
one as Ji„. Due to the linear dispersion relation near 
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the Fermi points, the densities are written by the sum 
of left and right current as Jout = JL,out + JR,out and 
Jin — Jlau + Jrau- Using these definitions of the left 
and right currents, we define a (quantum) field theoret- 
ical kinetic Hamiltonian respecting the linear dispersion 
as nil 
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where A — 2TThvp/L is the single-particle level spacing 
and is defined under the periodic boundary condition 
along a tubule axis. L is the length of the tube and 
vp is the Fermi velocity. The integral region is labeled 
as D G [0 : L]. 

The densities of each layer are expanded by the current 
operators as 
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The zero modes Qout and Qin denote the total charges 
in each layer respectively. The current operators satisfy 
the following bosonic commutation relations: 



And„ 



[J£, ( Jr)t] = AnSnr. 



(4) 
(5) 



and all of the other commutation relations vanish. Mak- 
ing use of the operator expression for the charge densi- 
ties, we rewrite the kinetic Hamiltonian as the sum of 
zero mode and non-zero modes. 



n>0 



(6) 



where 



A 

16 
A 

T 

A 

T 



)2 

' out 



^i.out) + Y5 ^^^^ ^ Ql^in) 1 C^) 



[(jr)tj£ + (j;j)tj«] 



(8) 



The operators Q^^out and Q^^m indicate the zero mode 
of the current in each layer, that is, they are equal to a 
integral of currents JL,out - Jr.oui and Jlau - Jrm over 
the nanotube length. 

The vacuum of this sector is defined as a direct prod- 
uct of the Dirac seas of both layers {\vac) — \vac,in) ® 
\vac,out)) which are labeled as \vac,in) and \vac,out) 
that satisfy the following conditions, 

J2\vac,in) = J^\vac,in) ~ 0, 
j2\vac, out) = j'^\vac, out) — 0, 



for positive n value. The vacuum state change into a 
new vacuum due to the intra and inter layer Coulomb 
interactions. They are given by 
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The Coulomb potentials for each layer are labeled as 
and y™. Between the layers, there is the Coulomb inter- 
action denoted V . Explicit form of these potentials would 
be difficult to observe experimentally. However we can 
drive some information of them from the stability of the 
systems, that we will discuss shortly. The potentials are 
rewritten by the Fourier series; 
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where the superscript a is a layer index and take an el- 
ement of the set {in, out, }. The Coulomb interactions 
are also decomposed into zero mode and non-zero modes 
as 



He 



(11) 



n>0 



where the non-zero modes are given by 
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[uiv+jr] [ji + UrV] 

Pn[{J2y + jR] [J2 + {JrV] 
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The Fourier components of the Coulomb interactions 
depend on system parameters; diameter and length of 
CNTs. Notice that /5™ and /3°"* for n ^ are modified 
by the vacuum polarization due to the massive bands, 
however the zero mode /J™ and do not receive any 
vacuum polarization correction. Standard Tomonaga- 
Luttinger theories of nanotubes take into account effects 
of non-linear interactions on low energy spectrum with 
an assumption (3^ = Pq [@, The Coulomb potentials 
and vacuum polarization depend on the wave number 
n. We respect the wave number dependent nature of 
the Coulomb interactions and vacuum polarization effect 
but do not consider the non-linear interactions including 
a tunneling interaction between layers. 



III. STABILITY CONDITION OF ZERO MODE 



The zero mode of the total Hamiltonian is 
Hp + H^ = ^— -I- -^o"*^ QoutQo 
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Suppose that the system is open, that is, number of the 
electron in each layer can be changed freely, then the total 
charges of each layer can take any value. Therefore when 
there is a direction in the space spanned by {Qout.Qin) 
that the energy of the zero mode become lowcir. the sys- 
tem is unstable along that direction. To avoid such in- 
stability, the energy must be elliptic as a function of the 
charges (Qout, Qin)- The condition of elliptic form corre- 
sponds to a negativity of the discriminant and results in 
an inequality Go < where 



Go = /3o - /3o" + 
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It is easy to find out that direction by diagonalizing the 

zero mode of the total Hamiltonian, to do so, we define 
new charges which are given by the following linear com- 
bination of the previous charges in each layer as 
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where A'' is a normalization constant. Making use of the 
new charges, we rewrite the zero mode of the Hamiltonian 
as 



-(- 

4 V 4 
1 /A 



+ 



Note that the condition of elliptic is equivalent to the pos- 
itively of the coefficient in front of Q^. Physical meaning 
of the new charges Q_ and Q+ is easily recognized in the 
case of /3o"* = In this case, we have 



(16) 



Hence, Q+ is a total charge mode and Q- corresponds to 
a capacitance charge. The system is unstable for a mov- 
ing electron between inner and outer layers when Go > 
and the system is open. 

On the other hand, when the system is closed, the total 
number of the charge is fixed, then the Hilbert space is 
restricted to satisfy the following condition: 



{Qin + Qout) \ closed) = 0. 



(17) 



The condition of stability of the system turns into a new 
condition, 



lA + i(/?g«* + /?j"-2/3o) 



>0. 



(18) 



It should be noted that the zero mode does not receive the 
vacuum polarization correction. Therefore, the condition 
for stability must be satisfied to keep the vacuum stable. 
However, the non zero modes are generally receiving the 
correction and depend on dispersion relations of massive 
bands. 



IV. STABILITY CONDITION OF NON-ZERO 
MODES 

We consider non-zero modes of the total Hamiltonian 
in this section. To analyze the energy spectrum of this 
system, first we apply the Bogoliubov transformations to 
the current operators of each layer separately and obtain 
the following expression of the Hamiltonian, 

+ an {Ji + {r^)^)x {01)^+31) 
+ «nOT + (iS)^)x((J2)UJS),(i9) 

where the transformed current operators are defined as 
follows, 

~jl \ _ ( coshC* sinht 
(j]^)tj - Vsinht,- 

J£ \ / coshij^" sinhi 

Prom the point of view of the new current operators, the 

interlayer Coulomb interaction is thought to have differ- 
ent coupling as compared with the original coupling /?„ 
that we denote a new interlayer coupling o„: 

a„ = /?„ (coshC* - sinhC*) (coshC - sinhC) . 

(20) 

The following Bogoliubov transformation angle and en- 
ergy spectrum were used to rewrite the non zero modes 
of the Hamiltonian, 

/3« 1 /A \ 

sinh2t« = cosh2f« = _ _ -f /3M , 



sinht°"* ^ 
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cosh C J 



i^" = ^yi + ^, a&{in,out}. 

Second we have to diagonalize the interacting terms of 
both currents whose coupling is a„. This can be done by 
shifting the current operator as 

11=11 + ^.{jl + iJ'li)'), (21) 

rH=rn+^Am+ji), (22) 

and the Hamiltonian leads to 

H^+H^ = K"* {01)12 + Qr)11) + ■ (23) 
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Here HI" is defined as 
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JTJout 
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This Hamiltonian implies that an effective attractive 
modes that couple to the Bogoliubov transformed cur- 
rent operators of the inner layer appears. The coupling 
depends on the energy spectrum of outer layer. Suppose 
that the outer layer is the semiconducting tube, then 
the excitation spectrum £'°"* is order of the gap energy 
which is very large as compared with the Coulomb en- 
ergy. Therefore the effective attractive interaction weak- 
ens much. 

We further diagonalize the term H™ as 

= E,, {{JiyJl + (J]^,)t jn) ^ (25) 
where new current operators are defined as 



J£ \ _ / cosh Sn sinh s„ \ / J£ 



{JrV ) coshs„y V (•^S)'^ 



(26) 



with the definition of the Bogoliubov transformation an- 
gle and energy spectrum 



sinh2s„ = — — , cosh2s„ — ^— El 
E„ ' V 
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From the energy eigenvalue of the lowest excitation 
mode, we read a condition of stability as 



En > 0, 



(27) 



which results in the negativity of the discriminant for the 
n-th Fourier component of the Coulomb interactions: 



A 



< 0. 



(28) 



This condition of stability is similar to that of the con- 
dition of stability of the zero mode for open systems. 
The physical entity of the lowest excitation mode is ca- 
pacitance mode between the layers. It is easy to prove 
this fact in the case of /?™ = In this case, to- 

tal density {p — Jin + J out) and capacitance density 
(tJ = Jin — Jout) decouple, so that one can calculate the 
energy spectrum of the capacitance modes without any 
difficulty. The energy spectrum is given by '^nE^ where 

E'^^ = ^\Jl + |-(/3r - Pn) and stability condition of the 
spectrum reduces to G„ < 0. 

We would like to consider several cases and check if 
the stability condition is satisfied. First, suppose that 
all the Coulomb potentials are equivalent, then /3°"* = 



Pti = gi'^ss G„ < which means the stability of 
the system. Second, assume the following conditions. 



/3n = 



/3r* + Pn 



oout Q}''^ 



then 



G„ = -|x ("1 + 2/3™) <0, 



(29) 



(30) 



which also result in the stability of the system. Finally, 
we analyze more realistic case. The diameter of the multi- 
walled nanotube is thick so that the vacuum polariza- 
tion due to the massive fermions loop can be no longer 
negligible and each Fourier component of the intra layer 
Coulomb potentials is modified. Therefore the condition 
of instability is given by 



Gn Pn 
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(32) 



We denote the static response functions as T™ and 
r°"* for each layer. Here we assume that the inter- 
layer Coulomb interaction does not receive the vacuum 
polarization corrections. This assumption is not correct 
for multi-walled CNTs with more than 3 metallic layers. 
When the vacuum polarization due to the massive bands 
is huge, the effective Coulomb couplings /3™ and /3°"* 
are close to zero. Therefore, the condition of instability 
reduces to the following condition: 



A 

Pn > TT- 



(33) 



We estimate the condition using the Fourier mode of the 
interlayer Coulomb potential which is given by 



Pn 



27rL7r 



21171(1 

L 



\ d 



(34) 



where Kq is the modified Bessel function, = 1.3[Ji] is 
a cutoff length which is introduced to take into account 
for the ionization energy of a 7r-electron and e is charge 
of electron which gives e^/An — 1.44[eV • nm]. d should 
be understood as the mean diameter of the outer and 
inner layer We plot the Fourier components of the 
Coulomb interaction and the energy A/8 in Fig.0. Note 
that there are some instability regions in the spectrum 
in this case. 

The above estimation in the case of vanishing intra- 
layer coupling is a rough estimation. Here we ana- 
lyze the stability condition using a static response func- 
tion derived at one-loop level P| supposing an equality 
/?™ = = (3n and that and received the 
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FIG. 1: Fourier components of the interlayer Coulomb in- 
teraction for several diameters. Horizontal axis denote wave 
number of the Coulomb potential where wave vector is de- 
fined by kn = 27rn/L. L is the length of the zigzag CNT that 
we take L = 3[^m] in these plots. 

same amount of the vacuum polarization correction. In 
this case, the condition of instabihty leads to 

Pn > /3r + V (35) 



We plot Pn and + ^ for several MWCNTs with dif- 
ferent diameter in Fig. ^. 




FIG. 2: Fourier components of the interlayer Coulomb inter- 
action and of the intra-layer interaction with one-loop vacuum 
polarization correction plus A/8. Horizontal axis denote wave 
number of the Coulomb potential, we take MWCNTs with 
L = 3[^m] and ds = 1.174, d_L = 7.045[nm]. Both metallic 
layers are assumed to be zigzag type CNTs. 

Figure. || shows that there are some unstable spectrum 
regions in momentum space, therefore we conclude that 



the system is thought to be unstable due to the vacuum 
polarization of the massive bands. Within the unstable 
region of wave number n Cz U, the excitation energy En 
(correctly 4n times En) take a negative value. Hence the 
system energy becomes lower if the current operators of 
the capacitance mode have a vacuum expectation value. 
Although we can not calculate the value of the amplitude, 
this instability would cause charge density oscillations in 
both layers. 



V. SUMMARY AND DISCUSSION 

In summary, the stability conditions on the Coulomb 
interactions in the multi-walled CNTs consisting of two 
metallic layers have been given within the framework of 
the simple model respecting the quantum nature of the 
electrons in linear bands and the Coulomb interactions. 
We have shown that because the massive bands screen 
the intra-layer Coulomb interactions in each layer, there 
arc finite unstable momentum regions in the system. Na- 
ture of the instability is governed by the capacitance 
mode between layers, therefore it is possible to detect 
the instability experimentally. 

The instability conditions are divided into zero and 
non-zero modes conditions. If the system is unstable, 
we can not predict a new ground state, but, can image 
several candidates of a new vacuum state. For example, 
if the zero mode is unstable, charge transfer between the 
layers is expected. After that charge transfer, the electric 
state of each layer looks like doped or un-doped single- 
wall CNTs (spontaneous capacitance). As a result of 
the moving of sufficient amount of charge between layers, 
the system may turn into a stable state because the low 
energy excitations are affected by the massive bands. On 
the other hand, even in the case of stable zero mode, 
the non-zero modes can still be unstable because of the 
vacuum polarization due to the massive bands. In this 
case, charge density oscillation in both layers may occur. 

We have not considered the tunneling interaction be- 
tween two layers in the present paper. Suppose that such 
an interaction is present, a change of the stability con- 
dition is expected. We would like to clarify this case in 
future work. 
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